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Abstract 

The microcanonical transfer matrix is used to study the distribution of the 
Fisher zeros of the Q > 2 Potts models in the complex temperature plane with 
nonzero external magnetic field Hq. Unlike the Ising model for Hq ^ which 
has only a non-physical critical point (the Fisher edge singularity), the Q > 2 
Potts models have physical critical points for Hq < as well as the Fisher edge 
singularities for Hq > 0. For Hq < the cross-over of the Fisher zeros of the 
Q-state Potts model into those of the {Q — l)-state Potts model is discussed, 
and the critical line of the three-state Potts ferromagnet is determined. For 
Hq > we investigate the edge singularity for finite lattices and compare 
our results with high-field, low-temperature series expansion of Enting. For 
3 < Q < 6 we find that the specific heat, magnetization, susceptibility, and 
the density of zeros diverge at the Fisher edge singularity with exponents ae, 
Pe, and 7e which satisfy the scaling law ag + 2/3e + 7e = 2. 
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I. INTRODUCTION 



The Q-state Potts model [jT],|^ in two dimensions exhibits a rich variety of critical behavior 
and is very fertile ground for the analytical and numerical investigation of first- and second- 
order phase transitions. With the exception of the Q = 2 Potts (Ising) model in the absence 
of an external magnetic field 0, exact solutions for arbitrary Q are not known. However, 
some exact results have been established for the Q-state Potts model. For Q = 2, 3 and 
4 there is a second-order phase transition, while for Q > 4 the transition is first order [Q. 
From the duality relation the critical temperature is known to be = J/kBln{l + a/^) [Q- 
For Q = 3 and 4 the critical exponents P| are known, while for Q > 4 the latent heat 0, 
spontaneous magnetization and correlation length at Tc are also known. 

By introducing the concept of the zeros of the partition function in the comp/ex magnetic- 
field plane (Yang-Lee zeros), Yang and Lee [§] proposed a mechanism for the occurrence of 
phase transitions in the thermodynamic limit and yielded a new insight into the unsolved 
problem of the Ising model in an arbitrary nonzero external magnetic field. Lee and Yang 
P] also formulated the celebrated circle theorem which states that the Yang-Lee zeros of 
the Ising ferromagnet lie on the unit circle xq = e*^ in the complex x = exp(/?if) plane 
for any size lattice and any type of boundary conditions. The density of zeros contains all 
the information about a system and in particular in the thermodynamic limit the density 
of zeros completely determine the critical behavior of the system H,^,^]. For example, the 
spontaneous magnetization of the Ising model is determined by the density of zeros on the 
positive real axis, i.e., mo(T) = 2'Kg{9 = 0,T). Above the critical temperature Tc, there is 
a gap in the distribution of zeros, centered at 6^ = 0, that is, g{0,T) = for \9\ < 6q{T). 
Within this gap the free energy is analytic and there is no phase transition. The Yang-Lee 
zeros at = i^o ^^re called the Yang-Lee edge zeros. As T — — > 0+, ^o(^) — * 0. At 
Tc the gap disappears, i.e., 6q{Tc) = 0, and g{0,Tc) = 0, which is the characteristic of a 
second-order phase transition. Below Tc, g{0,T) > and we have a finite spontaneous 
magnetization. Kortman and Griffiths Ijl^ carried out the first systematic investigation of 
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g{6,T), based on the high-field, high-temperature series expansion for the Ising model on 
the square lattice and the diamond lattice. They found that above T^, g{6,T) diverges at 
9 = ±6*0, i.e., at the Yang-Lee edge singularity for the square lattice. The divergence of 
the density of the Yang-Lee zeros means the magentization diverges, which does not occur 
at a physical critical point. Fisher |TT[] proposed the idea that the singularity at the Yang- 
Lee edge can be thought of as a new second-order phase transition with associated critical 
exponents and the Yang-Lee edge zero can be considered as a conventional critical point. 
Fisher also renamed the Yang-Lee edge zero as the Yang-Lee edge singularity. The critical 
point of the Yang-Lee edge singularity is associated with a 0^ theory, different from the 
usual critical point associated with the 0^ theory. The crossover dimension of the Yang-Lee 
edge singularity is = 6. The study of the Yang-Lee edge singularity has been extended to 
the classical n- vector model |jl2[, the quantum Heisenberg model [|1^], the spherical model 



T3| , the quantum one- dimensional transverse Ising model [jT4[, the hierarchical model [jT5[ 



and the one-dimensional Potts model |T6|. Using Fisher's idea and conformal field theory, 
Cardy |T^ studied the Yang-Lee edge singularity for a two-dimensional 0^ theory. Recently 



the Yang-Lee zeros of the two-dimensional Q-state Potts model have been studied ^1 



In 1964 Fisher |T^ emphasized that the partition function zeros in the complex tempera- 
ture plane (Fisher zeros) is also very useful in understanding phase transitions. In particular, 
in the complex temperature plane both the ferromagnetic phase and the antiferromagnetic 
phase can be considered at the same time. From the exact solutions ||^ of the square lattice 
Ising model Fisher conjectured that in the absence of an external magnetic field the zeros 
of the partition function lie on two circles in the complex y = exp{—f3J) plane given by 
Vfm = — 1 + V^e^^ (ferromagnetic circle) and uafm = 1 + V2e^^ (antiferromagnetic circle). 
Fisher also showed that the logarithmically infinite specific heat singularity of the Ising 



model results from the properties of the density of zeros. By numerical investigations |T9 
and analytical methods ||2^ it has been concluded that for very special boundary conditions 
the Fisher zeros of the Ising model do indeed lie on two circles, while for more general bound- 
ary conditions the zeros approach two circles as the size of lattices increases. Recently the 



locus of the Fisher zeros of the Q-state Potts model in the absence of an external magnetic 
field has been studied extensively pT|- p3| , p5| , pO| . It has been shown |^ that for self-dual 



boundary conditions near the ferromagnetic critical point r/c = + y/Q) the Fisher zeros 
of the Potts model on a finite square lattice lie on the circle with center —1/{Q — 1) and 
radius \/Q/{Q — 1) in the complex y-plane, while the antiferromagnetic circle of the Ising 
model completely disappears for Q > 2. It is also known that all the Fisher zeros of 



the one-state Potts model lie at y ^ = 0. Shrock et al. showed that for the two-dimensional 
Ising and Potts models in the absence of an external magnetic field there exist 



non-physical critical points in the complex temperature plane, at which thermodynamic 



functions including the magnetization diverge. Itzykson et al. considered the Fisher 
zeros in an external magnetic field for the first time. They studied the movement of the 
Fisher zero closest to the positive real axis for the Ising model as the strength of a magnetic 
field changes. For nonzero magnetic field there is a gap in the distribution of the Fisher 
zeros of the Ising model around the positive real axis even in the thermodynamic limit, 
which means that there is no phase transition. Matveev and Shrock studied the Fisher 
zeros of the two-dimensional Ising model in an external magnetic field using the high-field, 
low-temperature series expansion and the partition functions of finite-size systems. They 
found that for nonzero magnetic field the magnetization, susceptibility, specific heat, and 
the density of zeros diverge at the Fisher zero closest to the positive real axis, which we 
call the Fisher edge singularity. In this paper we discuss the Fisher zeros of the Q-state 
Potts model for nonzero magnetic field using the microcanonical transfer matrix and the 
high-field, low-temperature series expansion. 

II. MICROCANONICAL TRANSFER MATRIX AND SYMMETRIES 

The Q-state Potts model on a lattice G in an external magnetic field Hq is defined by 
the Hamiltonian 
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T<Q = -J - HgJ2S{ak,q), (1) 

<i,j> k 

where J is the couphng constant, < i,j > indicates a sum over nearest-neighbor pairs, 
CTj = 0, ...,Q — 1, and g is a fixed integer between and Q — 1. The partition function of the 
model is 

ZQ=Y.e-'''^, (2) 

where {cr„} denotes a sum over all possible configurations and (3 = {ksT)^^ . The partition 
function can be written as 

ZQ = y-'''T. Y.^QiE.M)x^y^ (3) 
E=a M=0 

where x = e^^'', y = e~^'^ , E and M are positive integers < E < Nf, and < M < Ng, 
respectively, Nfy and A^^ are the number of bonds and the number of sites on the lattice G, 
and Qq{E, M) is the number of states with fixed E and fixed M. Using the microcanonical 
transfer matrix (^uTM) p8|-31| we have calculated the number of states Qq{E,M) of the 



Q-state Potts model on finite square lattices with self-dual boundary conditions and 
cylindrical boundary conditions for 3 < Q < 8. 

In the absence of an external magnetic field the partition function of the Q-state Potts 
model is symmetric under the dual transformation 



iQ-i)y' 



(4) 



which gives the critical point 



and the invariant ferromagnetic circle of the Fisher zeros 



Voie) = • (6) 



The partition function of the Ising model has the additional symmetry 
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which maps the ferromagnetic Ising model to the antiferromagnetic model. This, together 
with the dual transformation, implies the invariance of the antiferromagnetic circle 

1/0 = 1 + \/2e^^ (8) 

However, the Q > 2 Potts models do not possess this second symmetry and the associated 
Fisher zeros are scattered in the non-critical region. For nonzero magnetic field the Ising 
model also has the symmetry 

x^-, (9) 

X 

and the Fisher zeros for < x < 1 have the same properties as those for 1 < x < oo. The 
Q > 2 Potts models do not have this symmetry, and distribution of zeros for < x < 1 
is different from the distribution for 1 < a; < oo. Because the Q > 2 Potts models are 
less symmetric than the Ising model, the zeros of the partition function have a much richer 
structure. For example, the Ising model has only non-physical critical points in the complex 
plane for x 7^ 1, while the Q > 2 Potts models have both non-physical and physical critical 
points in the same plane for x 7^ 1. In this paper we study the Fisher zeros of the Q-state 
Potts model for nonzero magnetic field to unveil some of the rich structures of the model. 

III. FISHER ZEROS OF THE THREE-STATE POTTS MODEL FOR x < 1 

In the limit Hq — > —00 {x — > 0) the partition function of the Q-state Potts model 
becomes 

Zq = y'"'' E m = 0)2/^' (10) 

E=Q 

where 0,q{E, M = 0) is the same as the number of states Qq_i(£') of the {Q — l)-state Potts 
model in the absence of an external magnetic field. As x decreases from 1 to 0, the Q-state 
Potts model is transformed into the [Q — l)-state Potts model in zero external field. Figure 
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1 shows the Fisher zeros in the complex ?/-plane of the three-state Potts model for x < 1 
with self-dual boundary conditions. Note that in the absence of an external magnetic field 
for self-dual boundary conditions the Fisher zeros in the critical region of the Potts model 
lie on the circle given by Eq. (6) In Figure 1 (a) the circle is that of Eq. (6) with 

Q = 3 (the three-state Potts circle), while in Figures 1 (c) and 1 (d) the circle is for Q = 2 
(the Ising circle). In Figure 1 (b) we show both the the three-state Potts circle (smaller one) 
and the Ising circle (larger one), and the Fisher zeros lie on neither the three-state Potts 
circle nor the Ising circle. In Figure 1 (c) the Fisher zeros near the ferromagnetic critical 
point begin to approach the Ising circle, and the antiferromagnetic circle of the Ising model 
begins to appear. In Figure 1 (d) almost all of the Fisher zeros, which will ultimately lie 
on the ferromagnetic circle of the Ising model at x = 0, are very close to this locus, and the 
antiferromagnetic circle becomes clearer. 

IV. CRITICAL POINT OF THE THREE-STATE POTTS MODEL IN A FIELD 

Hg<0 

For an external field Hg < 0, one of the Potts states is supressed relative to the others. 
The symmetry of the Hamiltonian is that of the {Q — l)-state Potts model in zero external 
field, so that we expect to see cross-over from the Q-state critical point to the {Q — l)-state 
critical point as —Hg is increased. 

We have studied the field dependence of the critical point for < x < 1 through the 
Fisher zero closest to the real axis, yi{x,L). For a given applied field yi approaches the 
critical point yi{L) — > ydx) in the limit L — > oo, and the thermal exponent yt{L) defined as 



22,26 



ln{Im[yi(L + l)]/Im[yi(L)]} 
= ln[(L + l)/L] ^''^ 

will approach the critical exponent yt{x). Table I shows values for ydx) extrapolated from 
calculations of yi{x,L) on L x L lattices for 3 < L < 8 using the Bulirsch-Stoer (BST) 



7 



algorithm |32|. The error estimates are twice the difference between the (n — 1,1) and 



(n — 1,2) approximants |^2|]. The critical points for x = 1 (three-state) and x = (two- 
state) Potts models are known exactly and are included in Table I for comparison. Note 
that the imaginary parts of ?/c(BST) are all consistent with zero. We have also calculated 
the thermal exponent, yt, applying the BST algorithm to the values given by Eq. (11), and 
these results are also presented in Table I. For x = 1 we find yt very close to the known value 
yt = 6/5 for the three-state model, but for x as large as 0.5 we obtain yt = 1, the value of 
the thermal exponent for the two-state (Ising) model. 

Figure 2 shows the critical line of the three-state Potts ferromagnet for Hg < 0. In Figure 
2 the upper line is the critical temperature of the two-state model, Tc{Q = 2) = l/ln(l-|-^/2), 
and the lower line is the critical temperature for the three-state model, Tc{Q = 3) = l/ln(l-|- 
V^)- The critical line for small —Hg is given by p3 



T ~ T,(Q = 3) ~ (-/7g)^'/^^ (12) 
where yt = 6/5 and yh = 28/15 for the three-state Potts model. 

V. FISHER ZEROS OF THE THREE-STATE POTTS MODEL FOR x>l 

In the limit Hg ^ oo {x —>■ oo) the positive field Hg favors the state q for every site and 
the Q-state Potts model is transformed into the one-state model ||23[. The zeros are given 
by 

Zq ~ y-""' E ^QiE, M = N,)y'' = 0. (13) 

_E=0 

Because VLq{E, M = Ng) = 1 for ii^ = and otherwise, Eq. (13) is 

= 0. (14) 

As X increases |?/o| for all the zeros increases without bound. 

Figure 3 shows the Fisher zeros in the complex y-plane of the three-state Potts model 
for h = (3H varying from to 4 in steps of 1. As /i increases, all the Fisher zeros move 
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away from the origin. Note that for h > there is accumulation of the Fisher zeros as we 
approach the Fisher edge zero, that is, the Fisher zero closest to the positive real axis. That 
kind of accumulation suggests that for h > the density of zeros diverges at the Fisher edge 
zero, which we call the Fisher edge singularity. The critical exponents associated with the 
edge singularity are defined in the usual way, 

(15) 

Ve 

me~(l-^)^% (16) 

Ve 

and 

Xe~(l-^)-^% (17) 
Ve 

where Ue is the location of the Fisher edge singularity, and Cg, nie and Xe are the singular 
parts of the specific heat, magnetization, and susceptibility, respectively. 

To study the critical behavior at the Fisher edge singularity we have used the high-field, 
low-temperature series expansion for the three-state Potts model due to Enting P^,|5B 



which is coded as partial generating functions. Table II shows estimates for t/e and (3e from 
Dlog Fade approximants [^] for the magnetization at a; = 100. For this value of x we find 
tte = 1.22(2), /3e = —0.197(6), and 7e = 1.20(3). Note that both /5e and are unphysical 
in that /^e < implies a divergent magnetization and ae > 1 implies a divergent energy 
density. The density of zeros near the Fisher edge singularity in the complex temperature 
plane is given by |jl8| 



(18) 

Ve 

Therefore, ae > 1 means that the density of zeros diverges at the Fisher edge singularity. 
From ae, Pe, and 7e we obtain 

ae + 2/?e + 7e = 2.03(4), (19) 
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so that the Rushbrooke scahng law a + 2/3 + 7 = 2, which is known to hold at a physical 
critical point, is also satisfied at the Fisher edge singularity. From the series expansions 
for the specific heat, magnetization, and susceptibility we have obtained the location of the 
Fisher edge singularity 

?/e(series) = 1.232(2) + 1.048(3)?, (20) 

which is in excellent agreement with the value we calculate by extrapolation from finite-size 
systems using the BST algorithm, 

ye(BST) = 1.233(7) + 1.050(3)z. (21) 

We have also studied the critical behavior at the Fisher edge singularity for several 
values of x. Table III shows the edge critical exponents and the location of the Fisher edge 
singularity for x = 20, 100 and 200. The edge critical exponents for any x satisfy the relation 

+ 2/5e + 7e = 2 withiu our error estimates. The locations of the Fisher edge singularity 
obtained from the series analysis agree very well with those extrapolated from finite size 
data by the BST algorithm. Table III suggests that the values of the edge critical exponents 
are independent of x. 

VI. FISHER ZEROS OF THE Q > 3 POTTS MODELS FOR NONZERO 

MAGNETIC FIELD 

Using the high-field, low-temperature series expansion of the Q-state Potts model for 
4 < Q < 6 we have studied the critical behavior at the Fisher edge singularity for 



Q > 3. Table IV shows the edge critical exponents and the locations of the Fisher edge 
singularities for 4 < Q < 6 and x = 100. The edge critical exponents for any Q satisfy the 
relation ae + 2j3e + 7e = 2 within our error estimates. As Q increases jSe appears to decrease 
slightly, while ae and 7e are constant within error. However, because the uncertainties in ae 
and 7e are large, we do not know whether ae and 7e are truly independent of Q. Even though 
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the Yang-Lee edge singularities have never been studied for the two-dimensional Q > 2 Potts 



models, according to the study of other models |TO|-p^,p^ and conformal field theory |[T] 
one expects the critical behavior of the Yang-Lee edge singularities in two dimensions to be 
universal. However, in a study of the Fisher (or complex-temperature) singularities of 
the Potts model in the absence of an external magnetic field Shrock et al. have observed 
a dependence of the edge critical exponents on Q. In Table IV the EST estimates and the 
series results for the location of the Fisher edge singularities agree with each other for Q = 4 
and 5. For Q = 6 we have calculated Qq{E, M) up to L = 5, and the EST extrapolation is 
unreliable because the maximum size of the lattice is small. Figure 4 shows the Fisher zeros 
in the complex y-plane of the six-state Potts model for x = 100, and the location of the 



edge singularities calculated from the series, which has been the traditional method p5|j27 
in the study of the Fisher (or complex-temperature) singularities. 



VII. CONCLUSION 

We have studied the Fisher zeros in the complex y-plane of the Q-state Potts model for 
X 7^ 1 using the microcanonical transfer matrix and the high-field, low-temperature series 
expansion. We have discussed the transformation of the Fisher zeros of the Q-state Potts 
model into those of the {Q — l)-state Potts model for x < 1, and into those of the one-state 
Potts model for x > 1. For a; < 1 we have obtained the critical line and calculated the critical 
exponents for several values of x. From the high-field, low temperature series expansion we 
have shown that for 3 < Q < 6 the specific heat, magnetization, susceptibility, and the 
density of zeros diverge algebraically at the Fisher edge singularity with characteristic edge 
exponents «£, Pe, and ■je- 
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TABLES 

TABLE I. The critical temperature yc and the critical exponent yt of the three-state Potts 
model for < X < 1. 



X yc (BST) yc (exact) yt (EST) yt (exact) 

0.414(3) + 0.0002(4)z 0.414213... 1.001(2) 1 
0.001 0.414(3) + 0.0002(4)i 1.001(2) 

0.05 0.413(5) + 0.0002(3)z 1.0009(6) 

0.5 0.400(2) + 0.000(2)z 0.982(21) 

1 0.366(2) + 0.0002(5)i 0.366025... 1.195(3) | 
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TABLE II. Values of /3e and ye 


estimated from Dlog 


; Pade approximants to magnetization for 


Q = 3 and x = 100. 






[N/D] 




Ve 


[7/8] 


-0.2077 


1.243256 + 1.052564? 


[8/8] 


-0.1967 


1.236604 + 1.045093i 


[9/8] 


-0.1952 


1.237445 + 1.043712i 


[8/9] 


-0.1954 


1.238006 + 1.043638? 


[9/9] 


-0.1957 


1.236733 + 1.044346? 


[10/9] 


-0.1939 


1.235093 + 1.044025? 


[9/10] 


-0.1956 


1.235556 + 1.044725? 


[10/10] 


-0.1953 


1.235416 + 1.044633? 


[11/10] 


-0.1930 


1.233087 + 1.045059? 


[10/11] 


-0.1956 


1.235523 + 1.044727? 


[11/11] 


-0.2003 


1.233911 + 1.047358? 


[12/11] 


-0.1975 


1.233556 + 1.046650? 


[11/12] 


-0.1972 


1.233675 + 1.046438? 


[12/12] 


-0.1977 


1.233538 + 1.046729? 


[13/12] 


-0.1975 


1.233549 + 1.046638? 


[12/13] 


-0.1959 


1.233478 + 1.046087? 


[13/13] 


-0.1991 


1.233249 + 1.047366? 


[14/13] 


-0.1966 


1.233340 + 1.046474? 


[13/14] 


-0.1943 


1.232810 + 1.046075? 


[14/14] 


-0.1937 


1.231850 + 1.046918? 


[15/14] 


-0.1970 


1.233458 + 1.046527? 


[14/15] 


-0.1945 


1.232889 + 1.046035? 


[15/15] 


-0.1945 


1.231888 + 1.047107? 


[16/15] 


-0.1966 


1.233328 + 1.046471? 
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TABLE III. The location and the edge critical exponents of the Fisher edge singularity for the 

three-state Potts model. 



X ye (BST) 


ye (series) 


ae 


Pe 


7e 


20 0.880(2) + 0.643(2)i 


0.880(4) +0.641(3)i 


1.2(1) 


-0.19(1) 


1.2(1) 


100 1.233(7) + 1.050(3)z 


1.232(2) + 1.048(3)z 


1.22(2) 


-0.197(6) 


1.20(3) 


200 1.436(8) + 1.270(l)i 


1.436(2) + 1.268(2)z 


1.22(2) 


-0.196(7) 


1.21(3) 


TABLE IV. The location and the edge critical exponents of the Fisher edge singularity for the 


4 < Q < 6 Potts models and x = 


100. 








Q ye (BST) 


ye (series) 


ae 




7e 


4 1.13(6) +0.96(4)z 


1.159(5) + 0.93(l)z 


1.18(8) 


-0.180(4) 


1.12(8) 


5 1.09(4) +0.861(8)i 


1.103(4) +0.86(l)i 


1.2(1) 


-0.173(4) 


1.1(1) 


6 


1.053(8) + 0.811(1)'/ 


1.2(1) 


-0.161(9) 


1.1(1) 
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FIG. 1. Fisher zeros in the complex y-plane of a 7 x 7 three-state Potts model for (a) x = 1, 
(b) 0.5, (c) 0.05, and (d) 0.001 with self-dual boundary conditions. 
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FIG. 2. Critical temperatures of the three-state Potts ferromagnet as a function of the magnetic 
field. Hq is in unit of J and T is in unit of J Iks- The upper dotted line is the Ising transition 
temperature in the limit Hq — > — oo, while the lower dotted line shows the critical temperature of 
the three-state Potts model for Ho = 0. 
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FIG. 3. Fisher zeros in the complex y-plane of an 8 x 8 three-state Potts model for h = f3H 
varying from to 4 in steps of 1 with cylindrical boundary conditions. 
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FIG. 4. Fisher zeros in the complex y-plane of a 5 x 5 six-state Potts model for x = 100 
with cylindrical boundary conditions. The two plus symbols show the locations of the Fisher edge 
singularities estimated from the series analysis. 



24 



